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Abstract 

Many large space system concepts will require active vibration control to satisfy crit- 
ical performance requirements such as line-of-sight pointing accuracy and constraints on 
root-mean-square (rms) surface roughness. In order for these concepts to become opera- 
tional, it is imperative that the benefits of active vibration control be practically demon- 
strated in ground-based experiments. This report describes the second experiment con- 
ducted by Harris as part of the NASA Controls-Structures Interaction (CSI) Guest Inves- 
tigator Program under contract NASl-18872. The results of this experiment demonstrate 
the successful application of the Maximum Entropy /Optimal Projection control design 
methodology to active vibration control for a flexible structure. The testbed is the Mini- 
MAST structure at NASA Langley Research Center and has features dynamically trace- 
able to future space systems. To maximize traceability to real flight systems the controllers 
were designed and implemented using sensors (four accelerometers and one rate gyro) that 
are actually mounted to the structure. Ground-mounted displacement sensors that could 
greatly ease the control design task were available but were used only for performance eval- 
uation. The use of the accelerometers increased the potential of destabilizing the system 
due to spillover effects and motivated the use of a precompensation strategy to achieve 
sufficient compensator roll-off. 
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1.0 INTRODUCTION 


Many future space missions will require active vibration control to satisfy critical performance 
requirements such as line-of-sight pointing accuracy and constraints on root-mean-square (rros) 
surface roughness. An important step in the development of this technology is demonstration 
in ground-based experiments. This' report discusses the second experiment conducted by Harris 
as part of the NASA Controls Structures Interaction (CSI) Guest Investigator Program. This 
experiment required control design and implementation for the Mini-MAST structure at NASA 
Langley Research Center and successfully demonstrated active structural control technology. The 
first experiment involved control design for the ACES structure at NASA Marshall Space Flight 
Center, Huntsville, AL and is described in [1-2]. 

The Mini-MAST structure, shown in Figure 1.1, is a beam-like truss structure. As demon- 
strated in Figure 1.2, this structure can be viewed as the secondary support tower of a precision 
optical structure. To achieve high accuracy line-of-sight pointing in the optical structure of Fig- 
ure 1.2 it is important to minimize the relative displacement of the tip of the beam with respect 
to the base. Hence, the primary objective of this experiment was to design controllers that pro- 
vide substantial reduction of the displacement of the tip of the Mini-MAST structure. Particular 
emphasis was also placed on controller simplicity (i.e., reduced-order and decentralized controller 
architectures). Complexity reduction in control law implementation is of paramount interest due 
to stringent limitations on throughput of even Btate-of-the-art space qualified processors. 

The primary methodology chosen for control design in this experiment was the Maximum En- 
tropy/Optimal Projection methodology [3-9], a subset of the Optimal Projection Approach for 
Uncertain Systems (OPUS) [10-11] which allows for the simultaneous trade-off of five fundamen- 
tal issues in control design: actuator sizing, sensor accuracy, controller order, robustness, and 
system performance. The Maximum Entropy/Optimal Projection approach was developed partic- 
ularly to enable the design of high performance, robust control laws for flexible structures. The 
design equations consist of four coupled matrix equations which specialize to the standard Linear- 
Quadratic-Gaussian (LQG) Riccati equations when the plant is known perfectly and a full order 
controller is desired. 

To maximize traceability to real flight systems, only the (acceleration and rate) sensors that 
are mounted on the Mini-MAST structure were used. Five sensors were used: four accelerometers 
and one rate gyro. Ground-mounted displacement sensors were available but were used only for 
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performance evaluation. Because of the two differentiators in the transfer functions from the 
control actuators to the accelerometers, the higher frequency modes were much more observable in 
the accelerometers than in the displacement sensors. Because the performance objective required 
control of the low frequency modes without destabilizing the higher frequency modes (a standard 
structural control problem), the use of accelerometers for control design significantly increased the 
spillover problem. Thus, in this case, it was much more challenging to achieve high performance 
design using accelerometers rather than displacement sensors. 

Because of the use of accelerometers, it was very important to ensure that the control laws 
rolled off sufficiently to avoid destabilizing the higher frequency modes. In this experiment the 
roll-off was enforced by using a precompensation strategy. That is, practical roll-off filters were 
first designed and included as part of the plant. The reduced-order LQG and Maximum Entropy 
control laws were designed using the modified plant. The roll-off filters were then appended to 
the reduced-order LQG and Maximum Entropy control laws to obtain the control laws which were 
actually implemented. As will be seen in the subsequent results, this methodology proved to be 
very effective for achieving the control design objectives for the Mini-MAST. 

The report is organized as follows. Section 2.0 provides a brief description of the Mini-MAST 
testbed while Section 3.0 describes the Mmi-MAST models provided by NASA Langley. Section 4.0 
reviews the Maximum Entropy /Optimal Projection approach to control design. Section 5.0 defines 
the basic control design problem and discusses the precompensation methodology used for control 
design. Section 6.0 describes the control design strategy used in this experiment and the reduced- 
order models that were used to develop the control laws. This section also discusses some of the 
key design details such as the dynamics of the precompensator filters and the weighting schemes 
that were used. Sections 7.0 and 8.0 describe respectively the decentralized and centralized designs 
that were implemented and the resultant performance improvement. Finally, Section 9.0 presents 
closing remarks and conclusions. 
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Figure 1.1 The Mini- MAST Test Article 
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Figure 1.2 The Traceability of the Mini-MAST Structure to Future Space Missions 



2.0 DESCRIPTION OF THE MINI-MAST TEST ARTICLE 


This section provides a brief description of the Mini-MAST experimental testbed, located at 
NASA Langley Research Center, Hampton, VA. A more detailed description is provided in [12]. 

The basic Mini-MAST test article is a generic space truss designed and manufactured by Astro 
Aerospace Corporation. The tubing members of the truss are made of graphite/epoxy. The truss 
beam is deployable and retractable and has a triangular cross section. The total height of the truss 
is 20.16 meters and the truss consists of 18 bays, each of which is 1.12 meters in height. 

The actuators and sensors available for control design implementation, disturbance generation, 
and performance evaluation are shown in Figure 1.1. The only actuators available for control are 
three torque wheel actuators that are mounted on the tip plate (top of Bay 18) parallel to the global 
x, y and z reference axes. The torque wheels provide both torsional and bending torque loads to 
the Mini-MAST. These DC permanent-magnet motors have a rated peak output of 50 ft-lbs at 50 
volts and 9.6 amps. 

The available control sensors are six Sundstrad Q A- 1400 servo accelerometers and three Watson 
angular rate gyros. Four accelerometers are located at the beam tip (Bay 18) and two are located 
on the mid platform (Bay 10). These sensors measure linear acceleration in the global x and y 
directions. The three rate sensors are located at the beam tip (Bay 18) and measure pitch (about 
the x-axis), roll (about the y-axis), and yaw (about the z-axis). 

Fifty-one Kaman KD-2300 proximity probes (i.e, displacement sensors) are installed on the 
support structure along the Mini-MAST. These devices can be used for control but were primarily 
intended for structural dynamic testing and performance evaluation. In our experimentation used 
the three Kaman sensors at Bay 18 for performance evaluation. The Kaman sensors at each bay 
have the orientation shown in Figure 2.1. 

Three Unholtz-Dickie 50-lb shakers are attached at Bay 9 for disturbance generation. These 
shakers are oriented normal to the faces of the truss at each of the three vertices as shown in Figure 
2 . 2 . 
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Figure 2.1 The Orientation of the Kaxnan Displacement Sensors 
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3.0 DESCRIPTION OF THE MINI-MAST FINITE ELEMENT MODELS 

Two models were provided by NASA Langley Research Center. The first model was used to 
generate the reduced order models that were used to design the decentralized controllers. A second 
model that had even better correspondence to experimental data was provided later in the program. 
This model was used to generate the reduced order models that were used to design the centralized 
controllers. The final evaluation model for each of the control designs was the full-order second 
model, discretized at 80 Hz, the sample frequency chosen for control law implementation. Delay 
states were included in the evaluation model to account for the computational delay. 

As shown in Figure 3.1, both of the models consisted of essentially three parts: the torque 
wheel dynamics, the finite element generated structural dynamics, and the 20 Hz analog Bessel 
filters used to filter the sensor outputs. The 20 Hz Bessel filters were chosen instead of lower 
frequency Bessel filters to allow better control of the second bending modes which have frequencies 
of around 6 Hz. The outputs and inputs for both models are listed in Table 3.2. Only Shaker-A was 
used to disturb the structure and Displacements-A, B and C of Bay 18 were used for performance 
evaluation. The inputs used for control law implementation were Torque-X, Torque-Y and Torque- 
Z while the outputs used were Acceleration- IX, Acceleration- 2Y, Acceleration-X, Acceleration-Y 
and Rate Gyro-Z. 

3.1 Model 1 

For this model the dynamics of the three torque wheels were given as follows. 


Torque Wheel X: 
Torque Wheel Y : 
Torque Wheel Z: 


Torque( Nm) 
Input Voltage 
Torque(iVm) 
Input Voltage 
Torque (Am) 
Input Voltage 


34227.1s 

(s + 22.8) (s + 407.6) 
38423.8s 

(s + 24.3) (s + 378.2) 
34293.0s 

(s + 24.6) (s + 370.5) 


The dynamics of the three torque wheels were quite similar to one another and significantly con- 


tributed to the plant dynamics over a wide bandwidth. For example, consider the Bode plots shown 


in Figure 3.3 of the dynamics of Torque Wheel X. It is seen that these dynamics provide 75° of 
phase lead at 0.8 Hz, the frequency of the first bending modes, and provide a magnitude variation 


of over 50 dB from 1 Hz to 100 Hz. 


The structural dynamics for this model were composed of 18 modes. The frequency, damping 
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and description of each mode is given in Table 3.5. The dynamics of each of the 20 Hz analog bessel 
filters were given as follows. 


20 Hz Bessel Filter: ^ U **P U * 

Input 


4.5475 x 10~ 18 (s a + 416.00s + 2.4275 x 10 18 ) 
(s + 208.23)(s* + 331.11a + 53012) 


NASA Langley was able to obtain good experimental data for the structure above .1 Hz and 
below 10 Hz. In this frequency regime Model 1 corresponded closely to experimental data. The most 
significant difference between this model and the actual system appeared to be due to inaccuracies 
in the dynamics of Torque Wheel X. The generally close correspondence is illustrated by Figures 
3.25 and 3.26 which were provided by NASA and compare the frequency responses of Model 1 
to the corresponding frequency responses derived directly from experimental data for two of the 
dominant transfer functions, Torque-X to Acceleration-2Y and Torque- Y to Acceleration- IX. 

3.2 Model 2 


The dynamics of the three torque wheels for this model differed somewhat from the dynamics 
provided as put of the first model and are given as follows. 


Torque Wheel X: 
Torque Wheel Y: 
Torque Wheel Z: 


Torque(iVm) _ 34861.3s 

Input Voltage ” (« + 23.50) (s + 336.90) 
Torque(Ym) _ 38508.9s 

Input Voltage (s + 23.56)(s + 401.31) 
Torque (Nm) _ 36433.0s 

Input Voltage (s + 23.44)(s + 372.34) 


The dynamics of Torque Wheel X differed the most from the corresponding dynamics of the first 
model. This difference is illustrated by the Bode plot comparison of Figure 3.3 and may be the 
most significant difference between the first and second models. The close correspondence between 
the current Torque Wheel X dynamics and experimental data is shown in Figure 3.4, which was 
provided by NASA. 


The structural dynamics for this model were composed of 28 modes. The frequency, damping 
and description of each mode is given in Table 3.6. The dynamics of the 20 Hz analog Bessel filters 
were identical to those given previously for the first model. 


Figures 3.7 thru 3.24 show the frequency responses of the transfer functions of Model 2 from 
each of the three torque wheel inputs to each of the five sensor outputs used for control law imple- 
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mentation, as weU as those from Shaker-A, used for disturbances, to the three tip displacements 
used for performance evaluation. Figures 3.8 and 3.12 also show that Model 2 is quite similar to 

Model 1. 

Model 2 corresponded closely to experimental data. The close correspondence is illustrated 
by Figures 3.27 and 3.28, which were provided by NASA and compare the frequency responses of 
Model 2 to the corresponding frequency responses derived directly from experimental data for two 
of the dominant transfer functions, Torque-X to Acceleration- 2Y and Torque-Y to Acceleration-lX. 
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Figure 3.1 The Three Parts of the Provided Models 







Table 3.2 The OutputB and Inputs of the Provided Models 
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Torque Wheel X Dynamics (Model 2 vs. Model 1) 






FREQ (HZ) 



Figure 3.4 Torque Wheel X Frequency Responses: Model 2 vs. Experimental 
(Model 2 = Solid Line) 
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Predicted 


Predicted 


Mode Number 


Frequency (Hz) 


1 

2 

3 

4 

5 

117 

118 

119 

120 
121 
122 

123 

124 

127 

128 

129 

130 

131 


.798 

.800 

4.37 

6.10 

6.16 

20.3 
21.6 

23.5 

28.6 
30.7 

32.1 

37.3 

38.3 

39.0 

42.2 
44.9 

54.3 

56.1 


Damping 

.018 

.018 

.012 

.010 

.010 

.005 

.005 

.005 

.005 

.005 

.005 

.005 

.005 

.005 

.005 

.005 

.005 

.005 


Description 

1st Y Bending 

1st X Bending 

1st Torsion 

2nd Y Bending 

2nd X Bending 

Tip Plate and Diagonals 

2nd Torsion 

Tip Plate, Bending and Diagonals 
Tip Plate, Bending and Diagonals 
3rd Y Bending 

3rd X Bending 

Mid Plate 

Mid Plate and Cable 

3rd Torsion 

4th Y Bending 
4th X Bending 

4th Torsion 

l ip Plate, Cable and Torsion 


Table 3.5 The Structural Modes of Model 1 


3-8 


















































Mode Number 
T7 
2 

3 

4 

5 

J 

7 

8 

9 

10 
33 

117 

118 

119 

120 
121 
122 
123 
126 

127 

128 
129 

132 

133 

134 

135 

136 

137 


Predicted 
Frequency (Hz) 


Predicted 

Damping 


Description 


0.83 


.018 


1st Y Bending 


0.83 


.018 


1st X Bending 


4.37 


.012 


1st Torsion 


6.38 


.010 


2nd Y Bending 


6.44 


.010 


2nd X Bending 


14.72 


.005 


Tip Plate 


14.83 


.005 


Diagonal 


15.42 


.005 


1st Axial 


15.57 


.005 


Tip Plate 


15.60 

17.17 

20.29 


21.80 


23.75 


.005 


.005 


.005 


.005 


.005 


Tip Plate 


Diagonal 


Tip Plate 


Second Torsion 


Both Plates 


25.83 


.005 


Both Plates 
3rd Y Bending 


31.62 


.005 


33.10 


.005 


3rd X Bending 


39.14 

42.46 


.005 


.005 


3rd Torsion 


4th Y Bending 


4th X Bending 


45.95 


.005 


54.14 


.005 


Tip Plate 


57.25 


60.85 

62.68 


.005 


.005 


.005 


4th Torsion 


Mid Plate 


Tip Plate 


66.04 


.005 


Tip Plate 


72.18 


72.37 


73.34 


.005 


.005 


.005 


5th Y Bending 


5th X Bending 


5th Torsion 


Table 3.6 The Structural Modes of Model 2 
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Torque X to Accel IX 



Torque X to Accel IX 



Figure 3.7 Bode Plots of Torque-X to Acceleration-lX for Model 2 
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Torque X to Accel 2Y 



Torque X to Accel 2Y 



Figure 3.8 Bode Plots of Torque-X to Acceleration-2Y for Models 1 and 2 
(Model 2 = Solid, Model 1 = Dashed) 
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Torque X to Accel Y 



Figure 3.10 Bode Plots of Torque-X to Acceleration-Y for Model 2 
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Torque Y to Accel 2Y 



Torque Y to Accel 2Y 



Figure 3.13 Bode Plots of Torque-Y to Acceleration-2Y for Model 2 
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Figure 3.11 
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Figure 3.18 Bode Plots of Torque-Z to Acceleration-2Y for Model 2 
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Figure 3.19 Bode Plots of Torque-Z to Acceleration-X for Model 2 
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Torque Z to Accel Y 



Figure 3.20 Bode Plots of Torque-Z to Acceleration-Y for Model 2 
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Figure 3.21 Bode Plot 


Torque Z to Rate Gyro Z 



Frequency (Hz) 
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of Torque- Z to Rate Gyro- Z for Model 2 
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Figure 3.22 
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Figure 3.25 Torque-X to Acceleration- 2 Y: Model 1 Bode Plots vs. Experimental 
Frequency Response Function (Experimental = Solid, 

Model 1 = Dashed) 
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Figure 3.27 Torque-X to Acceleration- 2Ys Model 2 Bode Plots vs. Experimental 
Frequency Response Function (Experimental = Solid, 

Model 2 = Dashed) 
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Figure 3.28 Torque- Y to Acceleration- IX: Model 2 Bode Plots vs. Experimental 
Frequency Response Function (Experimental = Solid, 

Model 2 = Dashed) 
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4.0 THE MAXIMUM ENTROPY/OPTIMAL PROJECTION APPROACH TO 
CONTROL DESIGN 

Maximum Entropy/Optimal Projection design [3—7] assumes an uncertain plant of the form 

x(t) = (A + AA)i(t) + (B + Afl)u(t) + u>i(t) (4.1) 

y(t) = (C + AC')i(t) + wj (t) (4.2) 


where * € IR"*, u e 01”“, y € IR"*, Wi € IR"* is white disturbance noise with nonnegative 
definite intensity Vi, tuj e IR"» is white observation noise with positive definite intensity Vj, and 
V)i and u >2 are uncorrelated. The state matrix A is assumed to be in real normal form such that A 
is block diagonal with scalar blocks corresponding to the real poles of the system and 2x2 blocks 

corresponding to the flexible modes. It is assumed that the uncertainty 


of the form 


-Vi W t 
-Wi -Vi 

A A is of the form 




A A = ^2 6*Ai> 1^*1 - Qt *> a » - 0 


(4.3) 


*=i 


where n a denotes the number of uncertain modes and each A,- is of the form 


Ai = block-diag{0, . . . ,0, 


0 1 

-1 0 


, 0 , . . . , 0 } 


(4.4) 


corresponding to uncertainty in the frequency of the i th mode. In this experiment, it was 
assumed that AS = 0 and A C = 0. 

Maximu m Entropy/Optimal Projection design allows the synthesis of a fixed order dynamic 
compensator, 


x c (t) = A c x c (t) + B e y(t) (4.5) 

u(i) = - C c x c (t ) (4.6) 

where x e G IR n ° and n c < n xy that (in some sense) minimizes the steady-state performance criterion 

J(A c ,B c ,C c ) = lim E[* T (t)Ri*(t) + u T S a u(0] (4.7) 

t— ►OO 

over the assumed uncertainty set. If we define 

E = BRi l B T and 2 = C T V 2 _1 C (4.8) 
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then the controller gains are given by 


A c = r(A, - E P- QE)G T 
B e = TQC^Vf 1 
C c = R^B t PG t 

where Q,P,Q and P are nonnegative definite solutions of the design equations 

0 = A.Q + QAj + V 1 -QtQ + f^ atAiQAf + £ a^QA? + r x QEQrJ 

1=1 »=1 

0 = Ajp + PAj + Ri - PEP + a? A? PA. + Y^a?A?PAi + rjFSFQr x 

»'=i *=i 

0 = (A, - EF)Q + Q(A. - E P) + QZQ - r x Q£QrJ 


0 = (A. - QE) t P + P(A, - QE) + PEP - r x PEPr x 

A A A A- - 

rank Q = rank P = rank QP = n c 
t = QP(QP)*, denotes the group generalized inverse 
and t has the factorization 


(4.9) 

(4.10) 

(4.11) 


(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 

(4.17) 


t = G t T, G,re nt n * xn -. ( 4 . 18 ) 

The Maximum Entropy design equations for full-order control law design are simply (4.12)-(4.15) 
with r x = 0. In this case r = G = / n> in (4.9)-(4.11). If we additionally assume perfect knowledge 
of the system (i.e., Ai = 0 for t = l,..,,n a ), then (4.12) and (4.13) become the standard LQG 
Riccati equations. 


The design equations (4.12)-(4.17) can be solved by using a homotopy algorithm [8,9]. As 
illustrated in Figure 4.1, this homotopy algorithm allows the deformation of an LQG controller into 
a full order Maximum Entropy controller. The Maximum Entropy controller is then reduced to an 
appropriate order by using an indirect controller reduction method. It is important that this initial 
reduced-order controller approximately solves the Maximum Entropy/Optimal Projection design 
equations to within a small error, although it is not even required to be a stabilizing controller. This 
can be achieved by beginning with a low authority LQG design and/or incorporating a sufficiently 
high level of uncertainty in the Maximum Entropy design. In practice a slight modification of the 
balanced controller reduction algorithm of Youssuff and Skelton [13] is currently used as the indirect 
controller reduction method. Once this initial reduced-order controller is obtained, the algorithm 
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is used to deform this controller into a Maximum Entropy/Optimal Projection controller. Then if 
a higher authority controller is desired, the final step of the algorithm is to increase the controller 
authority to a desirable level. 

In this experiment it was found that the final two steps were actually unnecessary. That 
is, for each design the weights of the initial LQG controller could always be chosen such that the 
initial reduced-order controllers had performance that could not be significantly improved by better 
satisfying the Maximum Entropy /Optimal Projection equations and were also of sufficiently high 
authority. This desirable feature of the design process was a nice by-product of the precompensation 
methodology discussed in the next section. 

Some key features of Maximum Entropy controllers are illustrated in Figures 4.2 thru 4.4 using 
controllers developed for the ACES structure. Each of these controllers used an approximately 
colocated angular rate sensor and torque actuator pair and were designed to attenuate all of the 
system modes less than 3 Hz. 

Figure 4.2 shows that in the performance region (DC-3Hz) Maximum Entropy design yielded 
positive real controllers providing robust stability to frequency uncertainties in the plant. Figure 4.3 
shows that, in the performance region, Maximum Entropy design smoothed out the compensator 
magnitudes, thus providing robust performance to shifts in the modal frequencies and indicating 
that the Maximum Entropy design reduced the order of the original LQG controller. 

Figure 4.4 shows that, outside of the controller bandwidth, the Maximum Entropy design ro- 
bustified the controller notches by increasing their width and depth. In the Mini-MAST experiment 
only this latter feature of Maximum Entropy design was needed since the LQG controllers were 
seen to be very robust within the controller bandwidth. This robustness was obtained primarily 
as a result of the wide separation between the first bending modes and the remaining structural 
modes and the precompensation methodology discussed in the next section. 
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’EASY PROBLEM' 



ti 


"EASY PROBLEM" 


Figure 4.1 The Practical Use of the Maximum Entropy /Optimal Projection 
Control Design Algorithm 
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PHASE IN DEGREES 


COMPENSATOR PHASE IN THE PERFORMANCE REGION 



Figure 4.2 An Illustration of a Phase Stabilization Property of Maximum 
Entropy Controllers 
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MAGNITUDE IN VOLTS/VOLT 



FREQ IN HZ 


Figure 4.3 An Illustration of Performance Robustness and an Order Reduction 
Property of Ma ximum Entropy Controllers 
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MAGNITUDE IN VOLTS/VOLT 


COMPENSATOR NOTCH 



FREQ IN HZ 


Figure 4.4 An Illustration of Ma ximum Entropy Notch Robustification 
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5.0 PROBLEM DEFINITION AND THE INTRODUCTION OF A PRECOMPEN- 
SATION METHODOLOGY FOR CONTROL DESIGN 

As mentioned in the Introduction, the basic control objective is to minimize the displacements 
at the tip of the Mini-Mast. It is assumed that the Mini-Mast is subjected to a .1 sec pulse of 50 N 
amplitude from Shaker-A. It is evident from Figures 3.19 through 3.21 that if Shaker-A is used to 
excite the Mini-Mast with a pulse of infinitesimal duration, then in the open-loop the displacements 
at the tip axe primarily due to the excitation of the first bending modes of the X and Y axes (at 
.8 Hz) with minor contributions from the first torsional mode (at 4.4 Hz) and the second bending 
modes of the X and Y axes (at 6.1 Hz and 6.2 Hz). Since the disturbance pulse is actually of finite 
duration, the contributions of the latter three modes is even further reduced. The dominance of 
the first bending modes in the open loop response is clearly seen in Figure 5.1 which shows the 
open-loop tip displacements to a .1 sec pulse from Shaker-A. 

From the above discussion it is apparent that the performance objectives require the control of 
no more than 5 modes with the first two bending modes being the two most important modes to 
attenuate. The control design thus requires the control of these 5 modes without destabilizing the 
remaining higher frequency modes, all of which are above 15 Hz. Recognize that the experimental 
data discussed in Section 3.0 that was used to verify the accuracy of the two models is only valid 
below 10 Hz. Thus, we do not have confidence in our models above 10 Hz and wish to ensure stability 
of these modes by having the compensators roll-off these modes to achieve gain stabilization. 

An additional motivation for using a roll-off strategy for achieving the desired high frequency 
stabilization is that we are actually going to use continuous-time models and design theory for the 
control designs. The controllers must then be discretized for implementation at 80 Hz. This dis- 
cretization will tend to distort the higher frequency dynamics of the continuous-time compensators 
which can cause instability when the discretized controllers are implemented. The use of compen- 
sator roll-off for high frequency stabilization reduces the importance of high frequency compensator 
dynamics. " 

It is important to note that, for this problem, the desired roll-off is not achieved by a straight 
forward application of LQG design theory. This fact is illustrated by Figure 5.2 which shows an LQG 
compensator designed using a reduced-order model of the Torque-X to Accel-2 Y transfer function. 
The reduced model included only one mode above 10 Hz, a 20.3 Hz mode. The LQG compensator 
was designed to attenuate the first bending mode at .8 Hz. The compensator shown in Figure 
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5.2 has sufficient gain at the first bending mode to achieve significant performance improvement 
but interacts with the 20 Hz mode by notching it and has enough gain above 20 Hz to potentially 
destabilize the unmodeled modes. This type of compensator is clearly undesirable. 

The desired roll-off can be ensured by assuming norm-bounded uncertainty and using either 
Hoo [13-18] or /i-synthesis [19,20] control design. Alternatively, one could use frequency weighted 
LQG [18,21,22] to attempt to achieve the desired roll-off. Here, however, we use a different strat- 
egy which, like frequency weighted LQG, is heuristic ally based. However, the “precompensation* 
methodology we used does have some attractive features. In particular, it avoids the great increase 
in dimensionality of the design model that is sometimes resultant from using Hoo, /i-synthesis or 
frequency weighted LQG methods. Also, this methodology Is very easy to implement. As illus- 
trated by the subsequent results, this method was very effective for achieving the control design 
objectives for the Mini-Mast structure. 

The basic idea behind the precompensation methodology is very simple and is illustrated in 
Figure 5.3 and 5.4. As shown in Figure 5.3, we simply embed the precompensation filters (in this 
case roll-off filters) , C u (a) and C y (a) in the plant a priori and design the Maximum Entropy /Optimal 
Projection controller £f(a) for this modified design plant. Then, as illustrated in Figure 5.4, the 
precompensation dynamics are included in the implemented compensator H(s). This methodology 
guarantees that the implemented controller will stabilize the nominal plant and also guarantees 
that the nominal performance will be that predicted in the design process since the closed loop 
eigenvalues of the feedback systems shown in Figure 5.3 and 5.4 are identical and for both figures 

— = + ^^22) ~ G12HG21 . . 

u; “ 1 + HG 22 

where 

H = C u HC y . (5.2) 

However, it is necessary to check after the design process as to whether the compensator achieved 
the desired roll-off (or robustness to high frequency unstructured uncertainty). 

The actual roll-off filters used are described in the next section. 
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Figure 5.3 Design Configuration for the Precompensation Methodology 
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H{s) = C u (s)H(s)C v (s) 


Figure 5.4 Actual Implementation for the Precompensation Methodology 
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6.0 CONTROL DESIGN STRATEGY AND MODELS 

A very practical control design and implementation strategy was used in this experiment. First, 
controllers with simple architectures (i.e., controllers which are decentralized, reduced order and 
use few sensors) were developed. Decentralized designs are more desirable for actual flight imple- 
mentation because they ease the processor requirements, but can be more difficult to design because 
decentralization requires a careful analysis of the physics of the structure. Subsequently, centralized 
control laws were developed in an attempt to improve the performance. The centralized design with 
the best performance did significantly improve the performance of the “best” decentralized design. 

From Figure 3.21 it was apparent that the transfer function from Torque-Z to Rate Gyro-Z 
was dominated by the first torsional mode at 4.4 Hz. Hence, it appeared feasible to use simple 
decentralized constant gain feedback from Rate Gyro-Z to Torque-Z to achieve high attenuation of 
the torsional mode. A full-order discrete-time model of the system was developed using Model 1 
to represent the system at the 80 Hz sampling frequency. This model was used to perform a root 
locus and it was determined from this root locus that the “optimal” gain was K=10. All of the 
implemented controllers were designed assuming that this feedback loop was closed. These designs 
added feedback loops that use only the accelerometers and the X and Y torques. The use of the 
Rate Gyro-Z to Torque-Z feedback essentially eliminated the influence of the torsional mode on 
the remaining loops as illustrated by Figure 6.1. The resultant performance improvement in the 
torsional motion is shown in Figure 6.2. 

Next, we considered attenuating the first bending modes in each axis by feeding back accelerom- 
eter outputs to the X and Y torque wheels. Since Acceleration-lX and Acceleration-2X provided 
essentially the same information about the acceleration of the Mini- MAST tip in the X-direction, 
we decided to use only one of these accelerometers and chose Acceleration- IX. Likewise we chose 
to use Acceleration-2Y instead of both Acceleration-2Y and Acceleration- 1Y. The remaining ac- 
celerometers, Acceleration-X and Acceleration-Y, were located near the center of the Mini-MAST 
at Bay 10. 

It was apparent from Figures 3.9, 3.10, 3.14 and 3.15 that the first bending modes were not 
very controllable and observable in the four transfer functions from Torque-X and Torque-Y to 
Acceleration-X and Acceleration-Y. However, Figures 3.7, 3.8, 3.12 and 3.13 revealed that the first 
han ding modes were very controllable and observable in the four transfer functions from Torque-X 
and Torque-Y to Acceleration- IX and Acceleration-2Y. Hence, we initially chose to feed back only 
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the latter two accelerometers to attenuate the first bending mode. Furthermore, it was evident from 
the latter four figures that the dominant transfer functions were from Torque-X to Acceleration- 
2Y and Torque- Y to Acceleration- IX. We thus chose to design feedback laws for the two loops 
corresponding to these transfer functions. 

Reduced-order continuous-time models were developed for the transfer functions Torque-X to 
Acceleration-2Y and Torque-Y to Acceleration-lX. First order all-pass filters were included in each 
model to account for the phase lag due to computational delay. Precompensation dynamics were 
added to each reduced-order model to yield the models actually used to design the Maximum 
Entropy /Optimal Projection controllers. The Torque-X to Acceleration-2 Y design model and the 
Torque-Y to Acceleration- IX design model each contained 20 states. 

The precompensation for each of the two control problems was chosen so that the modified 
design plant appeared to have a rate measurement and a true torque input within the control band- 
width. The precompensation thus consisted of dynamics to integrate the accelerometer output and 
dynamics to cancel the effects of the actuator dynamics below 10 Hz. This choice of precompensa- 
tion added the desired roll-off to the original design plant and also allowed us to develop a baseline 
rate- feedback design. Because the transfer functions for each of the two control problems were 
so similar, we decided to use the same precompensation for both. Referring to Figure 5.3, the 
precompensation used is given by 


C u {s) 


8 + 62.8 

a 2 + .628s + .3944 


Cy(B)= 1 . 


( 6 . 1 ) 


Due to the requirements of Maximum Entropy Design, each of the modified design models was 
placed in real normal form. Thus the A matrix of each of the modified design models was block- 


diagonal with one of the blocks of the form 


corresponding to the first bending mode. 


-v x 

-u /1 -v x 

The disturbance weighting matrix, Vi, and the performance state weighting matrix, R x> were each 


chosen to be block diagonal with the only nonzero block corresponding to the position of the first 

'0 0 


bending mode in the A matrix. The nonzero block of Vi is given by 


and the nonzero 


block of J?i is given by ^ where a is a scalar parameter which was used to determine the 
controller authority. This choice of weighting matrices was equivalent to assuming that the system 
was disturbed by a force actuator and the performance output was a displacement sensor that sees 
only the first bending mode when an impulse is placed in the system by the force actuator. Also, 
we chose i? 3 = Vj = 1. 


6-2 



The centralized designs used a reduced-order continuous-time model which contained 12 struc- 
tural modes. This model also included two first-order all-pass filters to account for the phase lag 
due to computational delay. Precompensation dynamics were added to the reduced order model 
which yielded a design model with 54 states. 

The precompensation dynamics that we used for the centralized controllers differed from that 
used for the decentralized controllers. The reason for this change was primarily an attempt to 
eliminate some of the low frequency oscillation that appeared in the output of the displacement 
sensors and ultimately limited the achievable performance of the decentralized designs. We initially 
attributed this low frequency oscillation to the migration of the low frequency poles of the precom- 
pensation dynamics towards the right half plane, although subsequent control designs and analyses 
revealed that this was not actually the case. Referring again to Figure 5.3, the precompensation 
used is given by 

C u (s) = diag{C Ut i(s),C Ul2 («)}, C y (s) = I (62) 


where 


39.0625 x 10 4 

C u ,i(«) = C Ui3 (e) - + 18 8755 + 6 2.5)(e J + 45.969 + 625) ' 


C U( i(a) and C u< j(e) are 4 Hz Butterworth filters. 


(6.3) 


The weighting scheme used for the centralized controllers is very straight-forward. The modified 
design model as before was placed in real normal form. Let u d £ IR denote the input corresponding 
to Shaker-A, z a € IR 3 denote a performance vector whose elements correspond to Displacements 
A, B and C of Bay 18, z h € IR 3 denote the performance vector whose elements correspond to 
Displacements A, B, and C of Bay 10, and let x denote the state vector of the modified design 
model in real normal form. We computed the input matrix D corresponding to u d and matrices 
Pa and Ph such that z a = P a x and z t = P b x. The disturbance matrix, V u and the state weighting 
were then chosen as 

V\ = aDD T , Rl = OtPa Pa + PPbPb ( 6 - 4 ) 


where the scalar parameter a was used to determine the controller authority. The scalar parameter 
P was chosen to be zero for the initial designs but was given a finite value in latter designs in 
order to penalize the Bay 10 displacements. The use of the parameter P was motivated by the fact 
that the dominant behaviour at Bay 10 was due to the second bending mode pair. Thus, P was 
essentially used to reduce the influence of the second bending modes on the displacement responses 
at Bay 18. 
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Figure 6.2 Open Loop vs. Closed Loop Rate Gyro-Z Response for Constant Gain 
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7.0 DECENTRALIZED CONTROL DESIGNS 

This section describes in more detail the decentralized controllers discussed in the previous 
section and presents the resultant performance improvement. Four controllers are described, each 
of which was designed using the precompensation methodology of Section 5.0. Each controller was 
also designed assuming that a decentralized constant gain feedback loop from Torque- Z to Rate 
Gyro-Z (with gain K=10) was closed. This feedback loop essentially eliminated the effects of the 
first torsional mode. Hence, although the controllers described used only Accelerometers- IX, 2Y, 
X and Y for sensing and Torques X and Y for actuation, it is important to keep in mind that 
they were implemented in conjunction with the constant gain feedback law from Torque-Z to Rate 

Gyro-Z. 

Each of the four decentralized controllers include at least two subcontrollers, one feeds back 
Accelerometer-2 Y to Torque-X while the other feeds back Accelerometer- IX to Torque- Y. These 
two subcontrollers were designed using the precompensation methodology of Section 6.0 with pre- 
compensation dynamics given by (6.1). Controller 1 consisted of two subcontrollers. Referring to 
Figure 5.3, these subcontrollers were designed by simply choosing H(s) equal to a constant. Since 
the precompensation was chosen to make the design plants C v (s)G(s)C u (s) for each subcontroller 
appear to have a rate measurement and a true torque input within the control bandwidth, Con- 
troller 1 approximates a rate feedback controller and was used as a baseline design. Controller 2 
also included only two subcontrollers and was designed without assuming uncertainty in any of 
the first five modes. Controller 3 was similar to Controller 2 except that it was designed assuming 
uncertainty in the second bending mode pair. Controller 4 consisted of Controller 3 plus two ad- 
ditional subcontrollers designed using classical control concepts. One of these subcontrollers feeds 
back Accelerometer- Y to Torque-X while the other subcontroller feeds back Accelerometer-X to 
Torque-Y. 

7.1 Controller 1: A Pseudo Rate Feedback Controller 

This controller consisted of two subcontrollers, one from Accelerometer-2 Y to Torque-X and 
one from Accelerometer- IX to Torque-Y. This controller was our least complex and, as previously 
mentioned, was used as a baseline design. The open loop vs. closed loop Bay 18 displacement 
responses are shown in Figure 7.1. Although this controUer did provide substantial increase in the 
performance, it will subsequently be seen that more complex controllers can achieve much better 
performance. 
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7.2 Controller 2; A Reduced Order LQG Controller Using Bay 18 Accelerometers 

Like Controller 1, this controller consisted of two subcontrollers, one from Accelerometer-2 Y 
to Torque-X and one from Accelerometer- IX to Torque- Y. However, the feedback laws in this case 
were designed using LQG. The decentralized feedback laws for each loop were each of order 9 
(including the precompensation) so that the implemented controller was of order 18. The open 
loop vs. closed loop Bay 18 displacement responses are shown in Figure 7.2. 

7.3 Controller 3: A Reduced Order Maximum Entropy Controller Using Bay 18 Ac- 
celerometers 

This controller was a robustified version of Controller 2 and was of order 20. It was designed 
assuming uncertainty in the second bending mode pair and, as shown in Figure 7.3, the Maxi- 
mum Entropy design robustified the corresponding notches for these two modes. This uncertainty 
description was originally motivated by an analysis of Controller 2 which indicated that the con- 
trollers were extremely sensitive to uncertainties in the second bending modes. However, due to 
the high fidelity of the models provided by NASA, subsequent analysis revealed that robustness 
was not really needed. Nevertheless, the Maximum Entropy design yielded useful gain margin at 
the frequency of the second bending modes and allowed us to scale the subcontrollers to further 
attenuate the influence of the second bending modes on the performance. The open loop vs. closed 
loop Bay 18 displacement responses are shown in Figure 7.4. 

7.4 Controller 4: Controller 3 plus a Classical Controller Using Bay 10 Accelerometers 

This controller consisted of Controller 3 plus two subcontrollers, H\(z) and H^(z) which were 
designed using classical control concepts. JIi(z) was used in the feedback path from Accelerometer- 
Y to Torque-X while Hi{z) was used in the feedback path from Accelerometer-X to Torque- Y. 
These two subcontrollers were designed to further attenuate the second bending modes and are 
given respectively by 

Hi(z) = K\zj[z — .9245) 

JJ 2 (z) = K^zjz{z — .9245) 2 

where K\ and Jfa are constant gains. The implemented compensator was of order 24. The open 
loop vs. closed loop Bay 18 displacement responses are shown in Figure 7.5. The corresponding 
displacement responses at Bay 10 are shown in Figure 7.6 while the commands to the torque wheels 
are shown in Figure 7.7. The desired attenuation of the second bending modes was achieved as can 
be seen by comparing the performance of Controller 3, as shown in Figure 7.4, to the performance 
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of Controller 4, shown in Figure 7.5. 

Figure 7.8 shows the performance for each of the controllers as evidenced from Displacement- 
A of Bay 18. The results were very close to the predicted results shown in Figure 7.9. This 
close correlation was largely due to the fidelity of the models provided by NASA Langley and is 
also partially due to the robustness of the controllers. Tha gains of Controller 4 are given in the 
Appendix. 
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BAY 18 (CONTROLLER 1) 



Figure 7.1 Open Loop vs. Closed Loop Bay 18 Displacement Responses for 
Controller 1 (Closed Loop = Solid, Open Loop = Dashed) 
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Figure 7.2 Open Loop vs. Closed Loop Bay 18 Displacement Responses for 
Controller 2 (Closed Loop = Solid, Open Loop = Dashed) 







Mag(dB) 



Frequency (Hz) 


Figure 7.3 Comparison of the Magnitude Bode Plots for the Acceleration^ to 
Torque-X Subcontroller of Controllers 2 and 3 (Controller 2 = Solid, 
Controller 3 = Dashed) 
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DISPLACEMENT B BAY 18 (CONTROLLER 3) 



Figure 7.4 Open Loop vs. Closed Loop Bay 18 Displacement Responses 
for Controller 3 (Closed Loop = Solid, Open Loop = Dashed) 








BAY 18 (CONTROLLER 4) 



Figure 7.5 Open Loop va. Closed Loop Bay 18 Displacement Responses 
for Controller 4 (Closed Loop = Solid, Open Loop = Dashed 
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Figure 7.6 Open Loop vs. Closed Loop Bay 10 Displacement Responses 
for Controller 4 (Closed Loop = Solid, Open Loop = Dashed) 




TORQUE X COMMANDS (CONTROLLER 4) 



Figure 7.7 Closed Loop Commands to the Torque Wheels for Controller 4 





The experimental and predicted open and closed loop responses of Displacement A of Bay 18 
for Controllers 1, 2, 3 and 4 are shown in the next two figures. 
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Figure 7.8 Open and Closed Loop Responses of Displacement A of Bay 18 
for Controllers 1 thru 4 


7-12 






Figure 7.9 
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8.0 centralize!) controllers 

This section further describes the centralized controllers discussed in Section 6 and presents 
the resultant performance improvement for each controller. Five controllers are described. As was 
the case for the decentralized controllers, Controllers 1 thru 4, each of these centralized controllers 
was designed assuming that a decentralized constant gain feedback loop from Torque-Z to Rate 
Gyro-Z (with gain K=10) was closed which effectively eliminated the effects of the first torsional 

mode. 

Each of the five controllers was designed using the precompensation methodology of Section 6 
with precompensation dynamics given by (6.2) and (6.3). The design weights are given by (6.4) 
with a and f) being the design parameters. The first two centralized controllers, Controllers 5 and 
6, feed back Accelerometers-2Y and IX to Torques-X and Y. Controller 6 is essentially a robustified 
version of Controller 5. Controllers 7 thru 9 also feedback Accelerometers-Y and X to Torques-X 
and Y. Neither of these latter three controllers was designed assuming uncertainty in the first five 
modes. The difference in these controllers was due to the amount of performance penalty that was 
placed on the displacements at Bay 10. 

It is of interest to note that although the centralized controllers were designed assuming 20 
Hz Bessel filters to process the sensor outputs, each of the controllers was originally inadvertently 
implemented with 10 Hz Bessel filters which provided substantial phase delay in the controller 
bandwidth. Nevertheless each controller was able to stabilize the system and obtain substantial 
performance improvement after the controller gain was reduced by a factor of 2. In addition, two 
of the controller, Controllers 7 and 8 were stabilizing without reducing their gain. 

8.1 Controller 5: A Reduced-Order LQG Controller Using Bay 18 Accelerometers 

This controller was the least complex centralized controller that we implemented and was 
designed to feed back Accelerometers-2Y and IX to Torques-X and Y. The displacements at Bay 
10 were not penalized in the design process (i.e., beta=0 in (6.4)). No uncertainty was assumed 
in the first five modes. The controller is of order 26 and the open loop vs. closed loop Bay 18 
displacement responses are shown in Figure 8.1. By comparing this figure with Figure 7.5, it is seen 
that the performance of this compensator, which uses two sensors, is comparable to the performance 
of Controller 4, a decentralized controller which uses 4 sensors and has better performance than 
the decentralized controllers using only the two sensors employed by Controller 5. 


8-1 


8.2 Controller 6: A Reduced-Order Maximum Entropy Controller Using Bay 18 Ac- 
celerometers 

This controller is a robustified version of Controller 5 and is also of order 26. It was designed 
assuming uncertainty in the second bending mode pair. The open loop vs. closed loop Bay 18 
displacement responses are shown in Figure 8.2. The performance was essentially the same as that 
of the performance of Controller 5, shown in Figure 8.1. 

8.3 Controller 7: A Reduced-Order LQG Controller Using Bay 10 and Bay 18 Ac- 
celerometers 

This controller was designed to feed back Accelerometers-2Y, IX, X and Y to Torques-X and Y. 
The displacements at Bay 10 were not penalized in the design process (i.e., /3 = 0 in (6.4)) and no 
uncertainty was assumed in the first five modes. The controller was of order 25. The open loop vs. 
closed loop Bay 18 displacement responses are shown in Figure 8.3. By comparing this figure with 
Figure 7.5, which shows the performance of Controller 4, it is seen that this centralized controller 
yielded significantly better performance than the highest performance obtained by a decentralized 
controller. 

8.4 Controller 8: Another Reduced-Order LQG Controller Using Bay 10 and Bay 18 
Accelerometers 

This controller used the same sensors and actuators as Controller 7 and also assumed no 
uncertainty in the first five modes. However, the controller did penalize the displacements at Bay 
10. In particular, /? and a in (6.4) were chosen such that /3/a = 2 in an attempt to further attenuate 
the second bending mode contribution to the displacement responses at Bay 18. The controller is of 
order 31 and the open loop vs. closed loop Bay 18 displacement responses are shown in Figure 8.4. 
By comparing this figure with Figure 8.3, it is seen that the closed-loop responses for Controller 8 
were slightly better than those seen for Controller 7. 

8.5 Controller 9: A Final Reduced-Order LQG Controller Using Bay 10 and Bay 18 
Accelerometers 

This controller used the same four sensors and two actuators as Controllers 7 and 8. The 
design process was essentially identical to that of Controller 8 except that a and j3 in (6.4) were 
chosen so that 3 / a = 5. The controller is of order 33 and the open loop vs. closed loop Bay 18 
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displacement responses are shown in Figure 8.5. The corresponding displacement responses at Bay 
10 are shown in Figure 8.6 while the commands to the torque wheels are shown in Figure 8.7. This 
controller yielded the best performance of any of the decentralized or centralized controllers that 

were implemented. 

The increasing performance improvement of Controller 4 (the highest performance decentral- 
ized controller) and Controllers 5, 7 and 9 is shown in Figure 8.8 by viewing the response of 
Displacement-A of Bay 18. The predicted performance of these three controllers is shown in Figure 
8.9. The close correspondence of these two figures is largely due to the high fidelity of the models 
provided by NASA Langley in the control bandwidth and is also partially due to the robustness of 
the implemented compensators. The gains of Controller 9 are given in the Appendix. 
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DISPLACEMENT B BAY 18 (CONTROLLER 5) 



Figure 8.1 Open Loop vs. Closed Loop Bay 18 Displacement Responses 
for Controller 5 (Closed Loop = Solid, Open Loop = Dashed) 





DISPLACEMENT B BAY 18 (CONTROLLER 6) 



Figure 8.2 Open Loop vs. Closed Loop Bay 18 Displacement Responses for 
Controller 8 (Closed Loop = Solid, Open Loop = Dashed) 









Figure 8.3 Open Loop vs. Closed Loop Bay 18 Displacement Responses 
for Controller 7 (Closed Loop = Solid, Open Loop = Dashed) 






DISPLACEMENT B BAY 18 (CONTROLLER 8) 



Figure 8.4 Open Loop vs. Closed Loop Bay 18 Displacement Responses 
for Controller 8 (Closed Loop = Solid, Open Loop = Dashed) 








BAY 18 (CONTROLLER 9) 



Figure 8.5 Open Loop vs. Closed Loop Bay 18 Displacement Responses 
for Controller 9 (Closed Loop = Solid, Open Loop = Dashed 










DISPLACEMENT B BAY 10 (CONTROLLER 9) 



Figure 8.6 Open Loop vs. Closed Loop Bay 10 Displacement Responses 
for Controller 9 (Closed Loop = Solid, Open Loop = Dashed) 
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Figure 8.7 Closed Loop Commands to the Torque Wheels for Controller 




The experimental and predicted open and closed loop responses of Displacement A of Bay 18 
for Controllers 4, 5, 7 and 9 are shown in the next two figures. 
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9.0 CLOSING REMARKS AND CONCLUSIONS 

This experiment successfully demonstrated high performance control law design and imple- 
mentation for the Mini-MAST testbed, a flexible structure which has features representative of 
future spacecraft. The controllers were designed using the Maximum Entropy /Optimal Projection 
approach to control design in conjunction with a precompensation methodology and thus provided 
validation of this strategy for control design for flexible structures. 

The precompensation methodology used in this experiment allows the designer to use classical 
insights to precompensate the original plant model at the inputs and outputs to develop a nice 
design plant for modern control design. The precompensation is appended to the modern controller 
to obtain the controller that is to be implemented. In this experiment the precompensation was 
chosen to provide high frequency roll-off (i.e., robustness with respect to unstructured uncertainty) 
and was also used to make the design plant appear to have rate sensors and torque inputs. This 
precompensation methodology allows the designer to obtain controller configurations that are very 
difficult or impossible to obtain using straight-forward applications of modern control theory. We 
believe that this strategy will be found to have other practical uses in the future and appears to 
be an effective means of integrating classical control concepts with modern control. 

The control design and implementation approach used was to start with simple controllers (i.e., 
reduced-order controllers, decentralized controllers and/or controllers using relatively few sensors 
and actuators) and increase controller complexity to increase performance. This strategy was ef- 
fective. The first four controllers each had a decentralized structure. The first controller was our 
simplest and did achieve some performance improvement. As the complexity of the controllers was 
increased by first adding order and then using additional sensors and controller order, the perfor- 
mance progressively improved. The performance was further improved by allowing the controllers 
to have a centralized structure. 

The results of this experiment and the Harris Guest Investigator experiment on the ACES 
structure [1,2] also illustrate that decentralized controllers and/or reduced-order controllers can 
provide very significant performance improvement for some flexible structure control problems. 
Controller simplicity is very important for the development of practical controllers due to substantial 
limits on throughput capability of even the most advanced space-qualified processors. 

The controllers in this experiment were designed to control the first five modes of the Mim- 
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MAST structure. The models provided by NASA Langley were quite accurate for these five modes 
and thus robustness within the controller bandwidth was not a real issue in this experiment as it was 
in the ACES experiment [1,2]. However, as already mentioned, the precompensation methodology 
was effective in providing the needed robustness for the higher frequency modes. 

The Harris team s experience in the NASA CSI Quest Investigator program and in imple- 
menting controllers on Harris testbeds [24] leads to the following recommendations for future flight 
experiments. 

1. Design the structure to allow decentralized control design (even if centralized controllers 
can yield better performance). 

2. Using the best model available without in-flight system identification, design highly robust 
decentralized controllers as the first controllers to be implemented. 

3. Design the mission to allow in-flight system identification using a proven identification tool 
such as the Eigensystem Realization Algorithm [25,26]. 

4. Using the system identification model, design a sequence of controllers of increasing com- 
plexity and potential performance improvement. 

5. Sequentially implement the sequence of controllers. If instability is sensed, then revert back 
to a controller that is known to be stabilizing. 

These steps should only be viewed as a rough outline. Many variations of this scheme are of course 
possible. 

Finally, one possible improvement is recommended in what is considered to be a very successful 
CSI Guest Investigator Program. It is recommended that future phases of this program be struc- 
tured to allow controllers to be implemented without requiring a guest investigator to be on-site. 
This would save the program significant cost and could thus allow more time for actual control 
design and analysis which would ultimately benefit NASA. 

Acknowledgement s. We thank Dr. Keith Belvin, Anne Bruner, Sharon Tanner, Jeff Sulla 
and David Geyer for their cooperation during this project. Without them this project would not 
have been successful. 
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11.0 APPENDIX CONTROL GAINS FOR CONTROLLERS 4 AND 9 


This Appendix presents the gains of Controllers 4 and 9 . It is assumed that each control law 

is of the form , . ... „ ... 

x e (k + 1 ) = A c x c (k) + B c y{k) 

u(k) = C c x 0 (k) + D c y{k). 

In addition the gains are presented in a basis in which the controller state matrix A e is in real 
normal form, i.e., 

A c = block-diag { A e , »'} £ 1 1 


where 


A c ,% = -Vi or A c< i = 



u>i 

~Vi 


and M < n e denotes the number of blocks. 

Controller 4 (A C 4,B C 4, C C 4,D c4 ) 
Description of Blocks of A c 4 


1. 0000000e+00 
2 . 0000000e+00 
3 . 00000006+00 
4 . 0000000e+00 
5 . 0000000 e +00 
6 . 00000006+00 
7 . 0000000e+00 
8 . 00000006+00 
9 . 00000006+00 
1 . 00000006+01 
1 . 10000006+01 
1. 2000000e+01 
1. 3000000 e +01 
1 . 40000006+01 
1 . 5000000 e +01 
1 . 60000006+01 
1.7000000e+01 
1 . 80000006+01 


size of A c ,i 

1 . 0000000 e +00 
2 . 00000006+00 
1 . 00000006+00 
1. 0000000 e +00 
1 . 0000000e+00 
1 . 00000006+00 
1. 0000000e+00 
2 . 00000006+00 
1 . 0000000 e +00 
1.0000000e+00 
2 . 0000000e+00 
2 . 0000000e+00 
2 . 00000006+00 
1.0000000e+00 
1 . 0000000 e +00 
1.0000000e+00 
2 . 0000000e+00 
1 . 0000000e+00 


Vi 

o.oooooooe+oo 

6 . 93639386-01 
8 . 11942236-01 
9 . 2446525 e -01 
9 . 24465256-01 
9 . 2446525 e -01 
9 . 3841864 e -01 
8 . 5729796 e -01 
9 . 6673736 e -01 
9 . 80687836-01 
8 . 7268108 e -01 
2 . 0521210 e -01 
2 . 0078435 e -01 
9 . 9605972 e -01 
9 . 9605972 e -01 
9 . 9605972 e -01 
9 . 9605972 e -01 
9 . 99933766-01 


Wi 

0 . 0000000 e +00 
7 . 1484760 e -02 
0 . 00000006+00 
0 . 0000000 e +00 
0 . 0000000 e +00 
0 . 0000000 e +00 
0 . 000000 0 e +00 
4 . 4416701 e -01 
0 . 00000 00 e +00 
0 . 00000006+00 
4 . 5172409 e -01 
9 . 7180085 e -01 
9 . 7505991 e -01 
6 . 7716165 e -03 
0 . 0000000 e +00 
6 . 7716165 e -03 

0.0000000e+00 
0 . 0000000 e +00 
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B C 4 = 


o. oooooooe+oo 
-4 . 783 1567 e+01 
-2 . 0131695e+02 
0.0000000e+00 
0 . 0000000e+00 
0 . 0000000e+00 
0 . 0000000e+00 
-2 . 4257444e+01 
-7 . 2516882e+00 
-4.3259987e+00 
0 . 0000000e+00 
-5 . 00028 53e+03 
0 . 0000000e+00 
0 • 0000000e+00 
-1.8654961e+00 
2 . 4169305e+00 
0.0000000e+00 
0 . 0000000e+00 
0.0000000e+00 
-6 . 5121393e+05 
0. 0000000e+00 
4 . 4661772e+06 
-2.0613049e+06 
0. OGOOOOOe+OO 


0.0000000e+00 
0 . 0000000e+00 
0.0000000e+00 
2 . 1982 125e+02 
0* 00G0000e+00 
0 . 0000000e+00 
0 . 0000000e+00 
0. 0000000e+00 
0. 00000Q0e+00 
0.0000000e+00 
1 . 2230791e+03 
0 . 0000000e+00 
-3 . 5172568e+00 
2 . 82084 69e+00 
0 * 0000000e+00 
0. 0000000e+00 
1.0533845e+00 
3 . 4778126e-01 
-1.9768639e+04 
0. 0000000e+00 
2 . 6873032e+06 
0. 0000000e+00 
0 . 0000000e+00 
-2 . 7080768e+06 


0. 0000000e+00 

o. oooooooe+oo 

0 . 0000000e+00 

0 . 0000000e+00 

0 . 0000000e+00 

0 . 0000000 e +00 

0 . 0000000 e +00 

0. 0000000e+00 

0.0000000e+00 

0. 0000000e+00 

O.OOOOOOOe+OO 

0. 0000000e+00 

0.0000000e+00 

0. 0000000e+00 

O.OOOOOOOe+OO 

0. 0000000e+00 

O.OOOOOOOe+OO 

O.OOOOOOOe+OO 

O.OOOOOOOe+OO 

O.OOOOOOOe+OO 

O.OOOOOOOe+OO 

O.OOOOOOOe+OO 

O.OOOOOOOe+OO 

O.OOOOOOOe+OO 



O.OOOOOOOe+OO 
9.2933761e-04 
2. 1761338e-03 
O.OOOOOOOe+OO 
9 . 2446525e-01 
O.OOOOOOOe+OO 
O.OOOOOOOe+OO 
-4.23229746-03 
4 . 1564 74 3e-03 
2 . 0288759e-03 
O.OOOOOOOe+OO 
-4 . 5837550e-03 
O.OOOOOOOe+OO 
O.OOOOOOOe+OO 
2.38492356-05 
6.2801479e-05 
O.OOOOOOOe+OO 
O.OOOOOOOe+OO 
O.OOOOOOOe+OO 
4 . 7453070e-03 
O.OOOOOOOe+OO 
-1.40526976-03 
-4 . 5329294e-03 
O.OOOOOOOe+OO 


1. 0000000e+00 
O.OOOOOOOe+OO 
O.OOOOOOOe+OO 
3.57160096-03 
O.OOOOOOOe+OO 
7 . 0710678e-01 
-7 . 0710678e-01 
O.OOOOOOOe+OO 
O.OOOOOOOe+OO 
0. 0000000e+00 
-4 . 963 608 Oe- 03 
O.OOOOOOOe+OO 
- 1 . 9571354 e -03 
9.9648177e-04 
o.oooooooe+oo 

O.OOOOOOOe+OO 
1 . 6183802e-04 
7 . 5593979e-05 
5 . 2 698732e-03 
0 . 0000000e+00 
-5 . 2730650e-03 
0. 0000000e+00 
0, 0000000e+00 
-5.2730669e-03 


0,0000000e+00 
0.0000000e+00 
0 . 0000000e+00 
0 . 0000000e+00 
0. 0000000e+00 
0.0000000e+00 
o. oooooooe+oo 
0 . 0000000e+00 
0 . OOOOOOOe+OO 
0 . OOOOOOOe+OO 
0.0000000e+00 
0.0000000e+00 
o.oooooooe+oo 

O.OOOOOOOe+OO 
0 • OOOOOOOe+OO 
O.OOOOOOOe+OO 
O.OOOOOOOe+OO 
0 . OOOOOOOe+OO 
0 . 0000000e+00 
0 . 0000000e+00 
O.OOOOOOOe+OO 
O.OOOOOOOe+OO 
0 . 0000000 e +00 
0. 0000000e+00 


0 . 0000000e+00 
0 . OOOOOOOe+OO 
0 . OOOOOOOe+OO 

o.oooooooe+oo 

2 • 0000000e-01 
0 • 00000 00e+00 
0 . 0000000e+00 
0 . OOOOOOOe+OO 
O.OOOOOOOe+OO 
0 . 0000000e+00 
O.OOOOOOOe+OO 
0 . OOOOOOOe+OO 
0 . 0000000e+00 
0. 0000000e+00 
0 • OOOOOOOe+OO 
0 • OOOOOOOe+OO 
0 . 0000000e+00 
0. 0000000e+00 
0. OOOOOOOe+OO 
0 • OOOOOOOe+OO 
0 . 0000000e+00 
0 • 0000000e+00 
0 . OOOOOOOe+OO 
0 • OOOOOOOe+OO 


Bc 4 = 


-4.07403566-02 
0 . OOOOOOOe+OO 
0 . OOOOOOOe+OO 


0. OOOOOOOe+OO 
2 . 5800693e-02 
0 . OOOOOOOe+OO 


O.OOOOOOOe+OO 
0 . OOOOOOOe+OO 
1. 0000000e+01 


O.OOOOOOOe+OO 
2 . 0000000e-01 
0 . 0000000e+00 


0 . 0000000e+00 
0. OOOOOOOe+OO 
0. OOOOOOOe+OO 
0. OOOOOOOe+OO 

o.oooooooe-*-oo 

1.86374106+15 
1 . 8637410e+15 
0. OOOOOOOe+OO 
O.OOOOOOOe+OO 
0.0000000e+00 
0 . 0000000e+00 
O.OOOOOOOe+OO 
0. OOOOOOOe+OO 
O.OOOOOOOe+OO 
O.OOOOOOOe+OO 
0 . OOOOOOOe+OO ; 

O.OOOOOOOe+OO 
0. 0000000e+00 
0. 0000000e+00 
0 . OOOOOOOe+OO 
0 . OOOOOOOe+OO 
0. OOOOOOOe+OO 
0. OOOOOOOe+OO I 
O.OOOOOOOe+OOJ 1 


O.OOOOOOOe+OO 
0 . 0000000e+00 
0 . 0000000e+00 
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Controller 9 (A C 9,B C 9 , C c q,D c q) 
Description of Blocks of A c 9 


l.oooooooe+oo 
2 . 0000000e+00 
3.0000000e+00 
4 . 0000000e+00 
5 . 0000000e+00 
6 . 0000000e+00 
7 . 0000000e+00 
8 . 0000000e+00 
9 . 0000000e+00 
1. 0000000e+01 
1. 1000000e+01 
1.2000000e+01 
1. 3000000e+01 
1 . 4000000e+01 
1 . 5000000e+01 
1.6000000e+01 
1 . 7000000e+01 
1 . 8000000e+01 


1 . 69299256+00 
-1.4306601e+00 
1 . 52281256+00 
. 3.2558851e+00 
-6 • 59387616-01 
-8.4178501e-02 
-2 . 1059136e+00 
-1 . 8 67 4 02 7 e+00 
-3 • 8143 185 e-01 
-9 . 169894 0e-01 
7 . 21619906-01 
-1 . 6G08374e+00 
6 . 0472042e-01 
3 * 2426651e-01 
-6 . 1423325e-01 
9 . 3971680e-02 
1 . 0151369e-03 
-1 . 9182462e-02 
8 . 1359307e-01 
1 . 9861650e-03, 
1. 3696630e+00 
3 . 0859232e-03’ 
-8 . 575 14 58e-03 
7 . 1495246e-03 
-3 . 2761330e-02 
-3 . 25007 88 e+01 
8 . 1864 104e+00 
-1 . 61663 30e-03 
-3 . 1116632e-04 
-1. 1474412e-03 
7 . 48053 02e-04 
-1. 3042720e-03 
3 . 7599438e-03 


size of A Ct i 

1.0000000e+00 
2 . 0000 000 e+00 
1.0000000e+00 
2 . 0000000e+00 
2 . 0000000e+00 
2 . 0000000e+00 
2 . 0000000e+00 
2 . 0000000e+00 
2 . 0000000e+00 
2 . 0000000e+00 
1. 0000000e+00 
2 . 0000000e+00 
2.0000000e+00 
2 . 0000000e+00 
2 . 0Q00000e+00 
2 . 0000000e+00 
2 . 0000000e+00 
2.0000000e+00 


2 . 0103 4 10e+00 
1.7241005e+00 
1 • 8122857e+00 
4 . 2949627e+00 
1 . 2435841e+00 
3.26053726-01 
4 . 3389398e+00 
1.8262086e+00 
6 . 8622161e-02 
2 . 2933 171e-01 
4 . 7232950e-01 
1.6817 193e+00 
6 . 4 15 1590e-02 
3.17326406-01 
3 . 7494903e-01 
1 . 9826487e-01 
7 . 2209392e-03 
1.31344046-02 
8 . 25638 10e-01 
7 . 4630388e-02 
9 . 8356352e-01 
4 . 3036453e-03 
1. 1573107e-02 
6 . 3036532e-03 
2 . 3846850e-02 
1. 5582181e+01 
4 . 0813536e+00 
6 . 8210711e-04 
2 . 0345945e-05 
7 . 288 1722e-04 
9. 1434180e-04 
4 . 5175594e-04 
■3 . 1484094e-03 


7 . 7466103e-05 
4 . 8671387e-02 
5 . 9485596e-01 
7.39837116-01 
7 . 39837 lle-01 
7.2259275e-01 
8.4674417e-01 
8 . 50482 17e-01 
8 . 50482 17e-01 
1 . 84 84 144e-01 
9 . 4688867e-01 
9 . 492 377 le-01 
-1. 8939909 e-02 
8 . 3618533e-01 
9 . 8927829e-01 
3 . 5073811e-01 
3 . 9195067e-01 
3 . 439563 3e-01 


5.1178601e+00 
-4.3819037e+00 
4 . 618 6803e+00 
-1 . 53 122 3 3 e+00 
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